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A detailed analysis of collective modes that couple to either vector or axial color currents in color- 
flavor locked phase of color superconducting dense quark matter at finite temperature is presented. 
Among the realm of collective modes, including the plasmons and the Nambu-Goldstone bosons, 
we also reveal the gapless Carlson-Goldman modes, resembling the scalar Nambu-Golstone bosons. 
These latter exist only in a close vicinity of the critical line. Their presence does not eliminate the 
Meissner effect, proving that the system remains in the color broken phase. The finite temperature 
properties of the plasmons and the Nambu-Goldstone bosons are also studied. In addition to the 
ordinary plasmon, we also reveal a "light" plasmon which has a narrow width and whose mass is of 
the order of the superconducting gap. 
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I. INTRODUCTION 

In the last few years much progress has been achieved in understanding the possible phases of quantum chromody- 
namics (QCD) at low temperatures and high densities. In particular it is expected that QCD at a sufficiently large 
density of quark matter possesses color superconducting phases . It is widely believed that one of color supercon- 
ducting phases may exist even at moderate densities that characterize the matter at the cores of compact stars (i.e., 
at densities just a few times larger than the density of the ordinary nuclear matter). The value of the superconducting 
order parameter is likely to be of order 10 to 100 MeV ||,^ (see also Refs. for estimates using the microscopic 

approach). These estimates of the order parameter are encouraging to seek for observable signatures from a color 
superconducting state inside compact stars While there are no many detailed studies in this direction yet, a 

couple of issues have already been addressed in Refs. p^-p^, and some of their conclusions seem to be very promising. 

At (asymptotically) large densities, the microscopic QCD becomes a weakly interacting theory This allows one 
to use analytical (although, non-perturbative) methods of quantum field theory for studying dense quark matter. At 
finite chemical potential /i, non-interacting quarks should occupy all the states with momenta \p\ < /i, thus, forming 
the Fermi surface. In a more realistic (weakly) interacting model, the attraction between quarks in the color-antitriplet 
channel leads to the famous Cooper instability. The latter is removed by a rearrangement of the ground state. As a 
result, in the case of the model with three quark flavors, the original gauge symmetry SU{3)c and the global chiral 
symmetry SU{3)l x SU{3)ii break down to the global diagonal SU{3)c+l+r subgroup ||2^. The corresponding phase 
is called the color-flavor locked (CFL) phase. Out of total sixteen candidates for Nambu-Goldstone (NG) bosons, 
eight are removed from the physical spectrum by the Higgs mechanism, providing masses to eight gluons. The other 
eight NG bosons show up as an octet (under the unbroken SU{3)c+l+r) of physical particles. In addition, the global 
baryon number symmetry as well as the approximate U{1)a symmetry also get broken. As a result, an extra NG 
boson and a pseudo-NG boson appear in the low energy spectrum. These latter particles are both singlets under 

SU{3)c+L+Fi.- rn m 

The low energy effective action of the CFL phase was derived in Refs. pi|-p5|, based on arguments of symmetry. 
Moreover, all parameters of the action were calculated in the limit of the asymptotically large chemical potential. 
The method of Refs. ^l|-p5|| is based on matching vacuum properties (such as vacuum energy and gluon screening) 
in the effective and microscopic theories. While being very powerful for many purposes, such a method is limited 
when it comes to determining the spectrum of bound states other than NG bosons. The other approach used to study 
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(diquark) bound states is based on the Bethe-Salpeter (BS) equation Ig^l- (A simplified study of the BS equation was 
also considered in Ref. and, for the case of two flavor dense QCD, in Ref. l2q].) 

In this paper, we study the spectrum of collective excitations in the CFL phase at finite temperature. To be more 
specific, we study only those collective modes that couple to the (vector and axial) color currents. The NG bosons 
and the massive plasmons are the examples of such modes. As we shall see, however, the plasmons and NG bosons do 
not exhaust all of the collective modes. In addition, we also predict the existence of the so-called Carlson-Goldman 
(CG) modes in the CFL phase of dense QCD (a short outline of this prediction was given in Ref. ||2^). In ordinary 
superconductors, the CG modes were experimentally discovered long time ago These modes appear only in a 

close vicinity of the critical line. Their quantum numbers are the same as those of the scalar NG bosons that are 
removed from the physical spectrum by the Higgs mechanism. Because of this, the appearance of the CG modes 
might look like a revival of the NG bosons |^ . 

This paper is organized as follows. In Sec. |ll[ we describe our model and introduce basic notations. Then, we 



present the general form of the effective action of the gluon field and discuss its most important properties in Sec. Ill 



In Sec. [V , we give the explicit form of the finite temperature improved hard-dense loop expression for the polarization 
tensor and reveal its internal structure. Then, by making use of the polarization tensor, we derive the general form of 
the current-current correlation function in Sec. M. In the same section, we derive the gluon and the (order parameter) 
phase field propagators and check that their poles appear at the same locations as in the current-current correlation 
function. In Sec. VI, we present a gen eral formalism for studying collective modes, coupled to vector (as well as 
In Sec. VII, VII] and IX, we analyze the properties of the plasmons, the NG bosons and the 



axial) color currents. 

gapless Carlson-Goldman modes, respectively. Sec. ^ deals with the Meissner efi^ect in the color superconducting 
phase. We derive the penetration depth of a constant magnetic field, and show that the superconducting quark 
matter corresponds to the so-called Pippard limit everywhere in the superconducting part of the phase diagram 
except for a very small region in the vicinity of the critical temperature. We summarize our results in Sec. XI. Some 
useful notations and formulas are collected in Appendices ^ and |^. In Appendix we present the general one- loop 
expression and the hard dense loop approximation of polarization tensor. The key integrals, used in the definition of 
the polarization tensor, are defined and approximately calculated in Appendices y and 



II. MODEL AND NOTATION 



As we mentioned in Introduction, the original SU{3)c x SU{3)l x SU{3)r symmetry of massless QCD breaks down 
to the global diagonal SU{3)c+l+r subgroup in the CFL ground state of dense quark matter. The color condensate 
in the CFL phase is given by the vacuum expectation value of the following diquark (antidiquark) field [ pO[ : 

(0| (^z,)^7' {^%)] |0) = ^iStS', + ^2S]Sl (1) 

where and = C"J'i3 are the Dirac spinor and its charge conjugate spinor, and C is a unitary matrix that 
satisfies C^^7^C = — 7J and C — —C^. In the last expression, we explicitly displayed the flavor = 1,2,3) and 
color (a,b — 1,2,3) indices of the spinor fields. The complex scalar quantities ki and K2 are determined by dynamics. 
Throughout the paper, we follow the notation of Ref. [g^ as close as possible. 

It is convenient to introduce the color-fiavor locked Weyl spinors ( octets and singlets under SU (3)c+l and SU{3)c+r, 
respectively) to replace the ordinary Dirac spinors. 











(2a) 










(2b) 










(2c) 
(2d) 



where A = 1, . . . , 8, and the sum over repeated indices is understood. Tilde denotes the charge conjugate spinors. 
Also, we use the conventional definition of the left- and right-handed projectors, 'P± = (1 ± 7^)/2. 

In the new notation, the non-zero order parameters ki and K2, as defined in Eq. (|^), are related to the following 
(singlet under the color-flavor locked residual symmetry) vacuum expectation values: 
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(OIV^^IO) = = ~l (3^1 + «2) , (3a) 

(oi^^v^^io) - -m^^^'m = ~ls^''^^2. (3b) 

In general, the order parameters ki and K2 are not independent. The details of the pairing dynamics should define 
a specific relation between the values of ki and K2- The value of their ratio, in its turn, defines a specific alignment 
of the order parameter. The analysis of the SD equation in phenomenological four-fermion models [Q, as well as 
in the microscopic QCD theory shows that the order parameter is dominated by the antitriplet-antitriplet 

contribution (ki — — K2)- While the sextet-sextet admixture (ki = K2) is small, it is never exactly zero. In fact, 
the latter does not break any additional symmetries |2^], and, for that reason alone, it cannot be excluded. In our 
analysis, we restrict ourselves to the case of a pure antitriplet-antitriplet order parameter. We do this simply for the 
purpose of convenience of presentation. The modification of our analysis to the case with a nonzero sextet-sextet 
order parameter is straightforward. In particular, it would only require to slightly change the ratio of the gaps for 
the quasiparticle in the octet and singlet channels (see below). 

The value of the order parameter in the CFL phase at zero temperature was estimated in Rcfs. j2^jlj,|l^, using 
phenomenological four-fermion as well as microscopic models [see Refs. [|| 0, dealing with the two flavor case]. While 
the exact value remains uncertain, it can be as large as 100 MeV. The temperature dependence of the diquark order 
parameter is essentially the same as that in the standard BCS theory of superconductivity JH]. Notice that the 
Meissner effect was neglected in the analysis of Ref. However, the recent result of Ref. |32f| suggests that such 

an approximation is justified. 

For completeness of presentation, we note that the inverse form of the relations in Eq. (H) reads 

i^n): = ^0^^ + + 7f + (4a) 

= ^{^^ + (^^)/ + 71 ('^ + (4b) 

Because the quarks in the color broken phase acquire a nonvanishing Majorana mass, it is also useful to introduce 
the left-handed and right-handed Majorana spinors, 

^^■ip + ip, *^ = ?A'^ + Vi^, (5a) 
$ = + 0, $-4 = 0^ + 0^, (5b) 

respectively. Then, in the CFL phase, the dynamically generated Majorana mass contribution takes a simple diagonal 
form, 

= (2AV- + 2AV+^ * - i*-^ (^AV- + AV+^ 

+ i$ (2AV+ + 2AV-^ $ - i|>^ {av+ + AV-^ (6) 
Here A = AyA+ -I- A^A~ and A — 7*'A^^7*', and the quark "on-shcU" projectors are 

" \P\ 



A± = - , 5 = 7"7. (7) 



Notice that the absolute value of the Majorana mass in Eq. (||) in the singlet channel is twice as large as the that in 
the octet channel. Also, the signs of those mass terms are opposite. Such a choice [see Eq. (^)] corresponds to the 
assumption of a pure antitriplet-antitriplet color-fiavor locked order parameter (ki = — K2)- 
By making use of the Majorana spinors in Eq. (||), the QCD action reads 

Cqcd = + M7°7') * + {i> + Ai7°7') + M7%') * + {i> " /n°7') 

4 2v6 
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Here the structure constants d^^'~' and f^^'^' are defined by the anticommutation and commutation relations of the 
Gell-Mann matrices, 



{A^,A^} =^<5^^ + 2rf^^^A^, (9) 

[A^,A^] =2z/^^c.A^. (10) 

By making use of the microscopic action in Eq. (||) , we could derive the well known Schwinger-Dyson equation for 
the quark propagator in the HDL improved, ladder approximation [l2| . Its solution reveals a dynamical generation of 
the Majorana quark masses. Such masses could be incorporated in the low energy quark action by the contribution 
Ca, presented in Eq. (||). 

In this paper, we are primarily interested in studying the properties of the collective modes at the non-perturbative 
vacuum in which the quarks have the Majorana masses. A specific mechanism of the mass generation is irrelevant 
for us here. Therefore, we could use the Hartree-Fock approach in our analysis. In particular, we add and subtract 
the mass term Ca to the microscopic action in Eq. (||). Then, we treat one of such terms as a part of the free quark 
propagator, while the other one as an interaction term. In this approach, a special care should be taken to preserve 
the gauge invariance of the model. This is because the Majorana mass term in Eq. (^) acquires phase factors under a 
general gauge transformation, and thus, it is not gauge invariant. To avoid the difficulty, it is convenient to consider 
the phases of the order parameter A as quantum fields restoring the gauge invariance of the model. 

By taking into account the structure of the dynamically generated mass terms in Eq. (^ , we arrive at the following 
quark propagators in the CFL phase, 

^ /(Po-e,)-2A^ j"(P0 + S)-2(A^)- 



for the left-handed fields. Similar expressions could also be written for the right-handed fields (the latter are obtained 
from the above expressions by the formal exchange of V- and V+). In the last two equations, we used the following 
notation: = |p| ± fji. As is clear from the structure of the quark propagators above, the value of the gap in the 
one-particle quark spectrum is determined by |A|p|. The so called antigap |A;J;|, on the other hand, has an obscure 
physical meaning because it can only be detected in the high energy (go ^ antiparticle excitations. In what follows 
we shall use the shorthand notation |At| = |A^| that should not create any confusion. 



III. EFFECTIVE ACTION OF GLUONS 



In this paper, we are mostly interested in collective modes, coupled to the (vector and axial) color currents. In 
the pursuit of a model independent treatment, we find it advantageous to start our consideration with the effective 
action, obtained by integrating out the quark degrees of freedom. Such an action could be easily derived in a wide 
class of four-fermion models of dense quark matter (in the weakly coupled regime) . As one could check, the effective 
Lagrangian density for the giuon field should have the following general structure: 

H^.(q)[<^-fig>f ]+..., (13) 

with the ellipsis denoting the interactions terms. Note that the presence of the phase field octet c/)^ in the last 
expression is very important for preserving the gauge invariance of the model. Under a gauge transformation, the 
phase in Eq. ( IT^ ) shifts so that it exactly compensates the transformation of the gluon field. Here we consider 
only infinitesimally small gauge transformations. In general, the gauge transformation of 0^ field is not a simple shift. 
However, it is always true that its transformation exactly compensates the transformation of the gluon field. 



'-'a,4> — d 



I?(°)(g) 



Ai 



A 
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The derivation of the effective Lagrangian density (y_3|) in dense QCD is not an easy task. In principle, one should 
start by integrating out the high energy gluons which, as is known from Refs. ||5|-[l^, play the dominant role in the 
dynamical generation of the gap. As a result of the integration, a new 4-fermion interaction term would be produced 
in the action. This latter could be treated by introducing the quantum Hubbard-Stratonovich composite field A 
with the quantum numbers of the color superconducting condensate. The effective Lagrangian density in Eq. (p^), 
then, is given by integrating out the quark degrees of freedom in the mean field approximation in which the vacuum 
expectation value of the composite field is fixed constant. This constant is identified with the value of the order 
parameter, (A) = A. The fields (j)^ are the phase fields of the composite field A that correspond to the (would be) 
NG bosons. 

While the microscopic theory of dense quark matter may truly be non-perturbative, the dynamics of gluon quasi- 
particles, given by the Lagrangian density in Eq. (|l^), could still be weakly interacting. This is known to be the case 
at least in the limit of asymptotically large quark densities. In this case, the explicit expression for the polarization 
tensor 11^;^ can be derived. Moreover, it is likely that the theory defined by Eq. is applicable even at realistic den- 
sities existing at the cores of compact stars (despite the fact that the derivation of Eq. (|l3|) fails at low density). This 
is because, in the case of weakly coupled gluons, the expression in Eq. (^) is the only possible quadratic expression 
consistent with the SU(3)c color gauge symmetry. 

In what follows, we use the simplest approximation for the polarization tensor fl^j, given by the improved hard 
dense loop (IHDL) approximation. Here "improved" indicates that the effect of the nonzero values of quark gaps are 
taken into account in the calculation of loop diagrams. The formal expression of such tensor was presented in Ref. . 
The zero temperature limit was also considered in Ref. |24| . Below, we present some other limits which has not been 
analyzed previously. In our notation, the explicit form of the (left-handed contribution of) IHDL polarization tensor 
(in the Matsubara formalism) reads 

+2Si{iujn,phfj.l'^Ss{iuJn + iVLm,p + q)lul^] , (14) 

where a;„ = (2n + l)7rr and r^„ i = 2m7r T ar e the fermion and boson Matsubara frequencies. In the derivation of 
Eq. (|lj), we used the identities (Al) and (A2) in Appendix^ The right-handed contribution H^.^^ is similar, and 



the full tensor is H^jy = H^^^^ -f H^^^j, (one could check, in fact, that I^R.^m = H^^^i^). 

Before going into further details, let us discuss the general properties of the polarization tensor in the IHDL 
approximation given in Eq. (p^. We start by emphasizing that such a tensor, obtained by integrating out the quark 
degrees of freedom, is not transverse in general (neither is it directly related to the spectral density of gluons, as we 
shall see later). This property is not connected with the non-Abelian nature of the system at hand. Indeed, to the 
leading order, dense QCD is essentially an Abelian theory, and many of its properties are similar to those in QED 
psf . Having said this, we should note right away that H^^ is transverse in the normal phase of the quark matter 
(above the critical temperature Tc) . This is because the NG phase fields are absent in the normal (symmetric) phase, 
and the polarization tensor is directly related to the observable current-current correlation function. 

In the broken phase (below Tc), on the other hand, a longitudinal part should necessarily appear. Of course, the 
existence of the longitudinal component in H^j^ cannot spoil the gauge invariance of the model in Eq. (^3|). To clarify 
this point, we should emphasize that the one-loop quark contribution LLi, is only one of two contributions to the 
complete expression for the gluon polarization tensor. As we show in Sec. ^ [see, for example, Eq. (^)], the other term 
comes from integrating out the "would be NG bosons" cfy^. A simple analysis shows that the final expression for the 
polarization tensor is in agreement with the Slavnov- Taylor identity required by gauge invariance. It is approapriate 
to mention that a similar observation at the level of the effective theory was also made in Ref. ISSl . 



As we show in Sec. VIII [sec also the arguments of Ref. [^], the longitudinal part of li^i, determines the dispersion 
relation of the pseudoscalar NG bosons in the broken phase. Some explanation is in order here. Strictly speaking, the 
properties of the NG bosons are related to the "axial" polarization tensor. In QCD at large density of quarks, this 
latter is approximately equal to the polarization tensor in Eq. (|l^). Thus, by using the leading order approximation, 
we could safely interchange the "axial" and "vector" polarization tensors. 

In passing, we note that the IHDL approximation, utilized in this paper, is selfconsistent. As one could check, 
the IHDL gluon propagator considerably deviates from the ordinary hard dense loop expression only in the infrared 
region go, |g| ^ A. This is the region where the Meissner effect is seen and where low energy modes come into the 
game. It is well known, however, that the relevant region of momenta in the gap equation is IS. <^ qQ\,q\ ^. Thus, 
using either the IHDL or ordinary hard dense loop approximation would not affect the leading order solution to the 
gap equation. At the same time, in order to preserve the gauge invariance of the model (e.g., in order to satisfy the 
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Slavnov- Taylor identity) it is crucial to use the IHDL approximation which, by construction, insures that the gauge 
invariance is intact. 



IV. GENERAL STRUCTURE OF THE ONE-LOOP POLARIZATION TENSOR 



By substituting the explicit form of quark propagators, given in Eqs. ( |lT| ) and (^h, into the definition of the 
polarization tensor in Eq. we derive the general result in Eq. (C2) in Appendix^Q Then, by dropping all the 

terms suppressed by powers of the chemical potential, we arrive at the following expression (here both the left- and 
right-handed contributions are taken into account): 



5V 



/ de 



u^u^il - ^^) + ^^g^^ + 



+ 



2\At\ 



2\An 



+ Upw^(i + r) .9m^ 



e^-^^l'?lV4-Wn(Wn+0m) 



e2|g]2/4-tj„K + n„) 



K + i?^)[K + ^mY + El] [ujI + i?2 )[(a,„ + r!„)2 + El] 



(c.2+ij2)[(^. 



-an)2+i?2]_ 

2iLj„e + iri„i(e - |g|C/2) 



|91 



(w2+i52)[K+a„)2+i?2] 



mm 



(1)^ 



(15) 



where we use the notation — (1,0,0,0) and ~ — (u ■ q)u^. Besides that, 0\iv{q) is one of the operators 



defined in Appendix see Eq. (Bl). This is the projector of the magnetic gluon modes. In Appendix we also 
introduce the projectors of the electric and unphysical (longitudinal in a 3-1-1 dimensional sense) modes, as well as 
intervening operator O^^uijq)^ mixing the electric and unphysical modes. It is appropriate to note here that the first 



term in Eq. (na) is the only term which got a contribution from the antiquark quasiparticles. The other terms came 



from low energy quasiparticles around the Fermi surface. These latter are characterized by the gap |At| 
the discussion at the end of Sec. ||). In Eq. (p^, we also use the notation: 



ig± = V >±^/2)^ + |AtP, 
^± = V(e±^l9l/2)2 + 4|ATP, 



(see 

(16) 
(17) 



which are the quasiparticle energies in the octet and singlet channels, respectively. In the calculation, we subtracted 
all divergencies independent of temperature and chemical potential. So, the result of Eq (15|) is already given in 
terms of renormalized quantities. Notice that the first term on the right hand side of Eq. (^) comes from divergent 
antiquark-quark loop contributions. In deriving the result in Eq. (|l5|), we performed the integration over the azimuthal 
angle, using the following relations: 



1 ^ 



1 - 3^2 q^q^ 

2 kT ■ 



(18) 
(19) 



As was already mentioned, in accordance with leading order approximation, we dropped all the terms suppressed by 
powers of the chemical potential in Eq. (p^). For example, we used 
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P+9 



(27r)3^'"'^ ^ (2 
\p + q]~fi^e- + + 0(|gf /m), 



(20) 
(21) 



where ^ is the cosine of the angle between p and q. Such an approximation is sufficient for studying the dynamics of 
the collective modes with \q\ <^ fi. Similarly, we approximated the results for the traces (with |g| <C |p| — fi) 



tr 



tr 
tr 



tr 



\p\' 



\p\ 



\P\ 



2 ■ 



(22) 
(23) 
(24) 
(25) 



The su m oy er t he M atsubara frequencies in Eq. ( [Tsl ) is straightforward. It reduces to three types of simple sums given 
in Eqs. (A7) - (A9) in Appendix By making use of the summation formulas also given in Appendix we derive 
another representation for the polarization tensor given in Eq. (C8). As is easy to check, then, the Lorentz structure 
of n^y(g) allows a well defined decomposition of the result in terms of the O^*^ operators, defined in Appendix 



n^.(g) = ni(q)OW(g) + n2(g)0(22(q) + n3((7)o(.t)(g) + n^{q)ol'J{q). 



(26) 



After a tedious, but straightforward rearrangement of different terms in the polarization tensor, we arrive at the 
following explicit results for the component functions: 



ni(g) - 




n2(<z) = 






n3(<z) = 


-I 


<il - W 


n4(<z) = 


-I 


ql - W 



[-\<fK{q) + qlL{q) + 2goklM(g)] , 

[qlK{q)-\cfL{q)~2qo\q\M{l)\. 
[qMK{q) - qML{q) - {ql + \q\^)M{q)\ 



(27a) 
(27b) 

(27c) 

(27d) 



where ujp — is the plasma frequency. The dimensionless scalar functions H{q), K{q), L{q) and M{q) on the 

right hand side have the following representation: 



H{q) 



12 



de 



4 



X{E-,E+)+ (2|ATp-7S-£;+-7e2 + 7-^^LL^ 



2\An 



E^E+ 



K{q) 



2\At\ + E^E+ - 



/ de 



4 



X{E^,E^ 
X{E^,E^ 



2|AtP - 7E-E+ + 7e' 



XiE-,E^ 



2\At\^ - E^E, 



^2 



X{E^,E^ 



2IA7 



E^E+ - 



2\At\^ - E^E+ -<? 



we 

4 

we 



YiE.,E+) 

Y{E^,E+) 
Y{E^,E+), (28) 



2|Ai 



7E-E+ + 7e^ -7 



\ee 



y{E-,E^ 



2IA7 



E^E^ 



^2 we 



Y{E^,E+) 
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2|AtP -E-E++e^ 



X{E-,E+)+ [2\ATr+E-E^ 



e^-^-^]Y{E-,E+), (29) 



2|AtP 



7E^E+ - 7e^ + | X{E^,E^ 



2|AtP - 7E^E+ - 



Y{E^,E^ 



2\At\^ + E^E^ 



2\At\'^ + E^E. 



e^ + \^]X{E^,E+) 



4 



2|AtP-S_£;+ -e^ + ^^L^ ) Y{E.,E+) 



X{E^,E^ 



M{q) - i J'^ J de 



7go 



E+-E- 
'•E,+E^ 



/ ,E^-E 



2 



X(i?_,i?+) + 7qo 



2|AtP -i?-£^4 



E+ + E- 



4 



E+ + E^ 



-90 



E+ + E^ 



XiE^,E+) + qo 





- E^ 


2 




E^ 




E+- 


E^ 


2 


E+ -f- 


E^ 




E+~ 


E_ 


2 



2 / 



where we introduced the foUowing two functions: 

(a + 6)[l-n(a)-n(&)] 

and n(a) and n(b) are the Fermi distribution functions. 



Y{a,b) 



Y{E^,E+) 
Y{E^,E+) 
Y{E^,E+ 

{a ~ b)[n{a) ~ n{b)] 
2ab[{a^b)^~ql] ' 



Y{E-,E^ 



(30) 



(31) 



(32) 



V. CURRENT-CURRENT CORRELATION FUNCTION AND PROPAGATORS 



The representation of the polarization tensor in Eq. ( |2^ ) in terms of three projectors and one intervening operator 
allows us to give a rather general treatment of the collective modes coupled to the vector and axial color currents. 
Before studying the properties of the collective modes in detail, it is instructive to consider the structure of the 
current-current correlation function in terms of the component functions of the polarization tensor Ili{q), see the 
definition in Eq. (p6|). 

Let us start with the Lagrangian density for the gauge field in Eq. (p^). In order to be able to generate the 
current-current correlation function, we introduce the classical source terms A^i'^{q) into the action, 

1 



A^'^iq) 



1 



Af^^{~q) + A^^^{-q) - iq'^q^^i-q)] n^.(g) U ^^"^(9) + A^''^ {q) + zg''0^(g) 



(33) 



where the ellipsis denote the interactions terms. Notice that the classical sources enter only through the second term 
of the Lagrangian density which is generated by integrating out the quarks. The gauge fixing parameter is denoted 
by A. As usual, it comes along with the projector of the unphysical gluon modes. The phase fields (f>^, required 
by gauge symmetry, enter the action in the same way as the longitudinal components of the gauge fields. Now, in 
order to derive the action for the transverse (physical) gluons which contain the information about the current-current 
correlation functions, we should integrate out the phase fields (j)^. Thus, we obtain the following expression: 



--A^-n-g)<Z^ [0W(g)+0(2j(^) + A0(^)(g)J A^''^{q) 



A^^^(-g)+A^'^(-g) 



n^.(9)- 



qpUP'^qa 



Af,'''{q) + A-'-^{q) 



(34) 
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By substituting the decomposition of the polarization tensor n^i/ in Eq. (M) , we arrive at the action: 



'A^^^{-q)+A^'^^{-q) 



ni(g)oJ,l'('z) + (n2(g) 



mq)\' 



0^{q) 



A^^''{q)+A^'%q) 



n3(g) 

Finally, by integrating out the quantum gauge fields, we arrive at the following quadratic form: 



i-q) 



q^n 



q2 [U2U3 + {u,y 



Hi 



{q' 



rO^'Mq) 



A^r 



(35) 



(36) 



As is clear, the tensor in the square brackets of this quadratic form defines the current-current correlation function, 
i.e.. 



n 



q^ [n^ns + [w^f] 
(g2 + n2)n3 + (n4)2 



0^{q) 



The location of its poles. 



q^ + Hi [q) = 0, "magnetic" , 
[q^ + n2 {q)] Ha {q) + "electric" , 



(37) 



(38) 
(39) 



defines (in a gauge invariant way) the spectrum of collective excitations in the model at hand. It is interesting to 
point out that one could have extracted this spectrum of the collective modes already from the poles of the quantum 
gauge field propagator in Eq. (|3^): 



q^ 



-^o«(,) 



^^3 



(q2 + n2)n3 



-(n^O(l)(,)-3^0(3)(,). 



(40) 



We performed the last integration just to show explicitly that the proper definition of the current-current correlation 
function indeed gives the anticipated result. 

Before concluding this section, it is also instructive to perform the following exercise. Let us perform the integration 
in Eq. (§3|) over the quantum gauge field first. Then, the intermediate action for the "phase" field reads: 



'^^''^(-?)-*9^'r^(-?) 
\q^ [(g2+ 02)03 



q^n 



Iromq) 



g2 [(Ag2+ 03)02 



(n4)2] 



Hi 



Ag2(g2+ 02) + [(g2+ 02)03 -f(n4)2] 



Kq^fiii 



Ag2(q2+02) + [(g2 + 02)03 + (n4)2] 



0^{q) 



Ag2(g2 + 02) + [(g2 + 02)03 + (04)2] 



At^\q)+iq^n-{q) 



0^}{q) 



(41) 



At first glance, it does not make sense to put any physical relevance behind the propagator of the field, because 
the latter is a gauge dependent phase field. In condensed matter physics, however, this type of propagator is often 
used to extract the spectrum of the electric-type collective modes, such as the CG modes. The explicit form of the 
(j)^ field propagator reads 



02) + [(g2 + 02)03 + (04)2] 



A(<72)2[(q2+n2)03 + (04)2] 



(42) 



We see that the location of one of the poles gives exactly the same spectrum as that in Eq. ( |39| ) . It is remarkable that 
the position of the corresponding pole and its residue are gauge invariant. Oowever, there exists another pole in the 
propagator in Eq. (^) at — Q. This latter is unphysical in a sense that it does not appear in the current-current 
correlation function, given in Eq. (|36|). The appearance of the unphysical pole clearly shows that the use of the 
propagator (42) for extracting spectra of collecting modes is limited. For completeness, let us also mention that there 
is no trace of the magnetic collective modes in the propagator of the longitudinal cj)"^ field. The information about the 
magnetic modes, however, is not lost in Eq. (p|). As is easy to check, after integrating out the (j)^ field in Eq. (|4]]), 
we arrive at the same result as in Eq. (^). 
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VI. COLLECTIVE MODES: GENERAL FORMALISM 



In the preceding section, we derived the general expression for the current-current correlation function, see Eq. (|37 
The location of its poles gives the spectra of the electric and magnetic type collective modes (assuming that their 
width is not too laxge). The explicit form of such spectra is determined by Eqs. ( ^8|) and (^). By making use of the 
definition in Eq. (p7|), the electric and magnetic dispersion relations take the following form: 



ql - + ^lH{q) = 0, "magnetic" , 
qlK{q) - \q\^L{q) - 2q^\q\M{q) + c^^ [K{q)L{q) + M^q)] = 0, "electric". 

Here we used the following identity, relating different component functions: 



n,n 



(n4 



ol [KL + Af 2] 



(43) 
(44) 

(45) 



To the leading order in powers of go E^nd \q\ (keeping the ratio v = qo/\q\ arbitrary and |At| finite), we obtain the 
following result for the component functions of the polarization tensor: 



F(g) = -l + i / 



9|AtP tanh(/3£;/2) (5|AtP + l^e^)n' {E) 



2E^ 



1 - 



3E\I SEqp 



1 - 



Eqo 



(2|AtP -e^ + EE){E + E)[l - n{E) - n{E)] (2|AtP - - EE){E - E)[n{E) - n{E)] 



K{q) 



de 



EE{E + EY 

5|ATptanh(/3S/2) 
4B3 ^ 



lAe^)n'{E) 



2E'^ 



EE{E~EY 

Eqo Eqo + e\q\ 
2e\q] "^Eqo-elq] 



(46) 



(2|AtP + - EE){E + E)[l - n{E) - n{E)] (2|AtP + + EE){E - E)[n{E) - n{E)] 



de 



EE{E + EY 
9|AtP tanh(/3£:/2) (5|AtP + Ue^)n'{E) 



AE^ 



2E'^ 



(2|AtP - + EE){E + E)[\ - n {E) - n{E)] ^ (2|At|^ - e 



EE{E-EY 

^E^ql iE\l Eqo + e\q\ 
e^W 2e3|g13 ^^Eqo-e\q\ 

2 EE) {E - E) [n{E) - n{E)] 



(47) 



M{q) 



7qo_ 
3 191 



EE{E + EY 
2e\q\ Eqo - e\q\ 



EE{E-EY 



n'{E), 



(48) 
(49) 



where the quasiparticle energies in the octet and singlet channels are E = V^^ + I^Tp and E ^ y?T4|A^, 
respectively, and (3 = 1/T is the inverse temperature. 

For further convenience, we rewrite these representations in the following shorthand form: 



H{v) 
K{v) 
L{v) 
M{v) 



-1+g [9/1+/3 



3 [-5/l+/2+/4 + /6(v)] 



-1+g [9/1+ /3 

|/s(.), 



/7(^')], 



(50a) 
(50b) 
(50c) 
(50d) 



where the li functions on the right hand side are defined in Appendix |^ In the same Appendix, we also present 
their nearcritical asymptotes which include all corrections up to the second power in |Ar|/T. 



VII. PLASMONS 



Dense quark matter is a non-Abelian plasma. Similarly to the ordinary plasma, it is characterized by the so- 
called plasma frequency tOp. By definition, the properties of the gluon field excitations change qualitatively when the 
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frequency increases from the values less than ujp to the values greater than LUp. For example, in ordinary plasma (in 
absence of magnetic fields) , only the modes with frequencies greater than the plasma frequency could freely propagate 
in the bulk. Such modes are called plasmons. 

To derive the spectrum of the plasmon excitations, it is sufficient to consider only the long wave length limit |gl — > 0, 
keeping go large. The only restriction is that qq <^ fi. As we shall see, such an assumption is justified in the weakly 
coupled limit. 

Let us start our analysis by writing down the component functions, defined in Eqs. ( p8| ) - (|3l|), in the long wave 
length limit. By substituting \q\ = 0, we obtain 



H{qo)=L{qo) = -l + ^ I de 



9|AtP tanh(/3£;/2) (2|AtP - + EE){E + E)[l - n{E) - n{E)] 



E{AE^ - ql) 



EE[{E + EY-ql] 



(2|A, 



- EE){E - E)[n{E) - n{E)] 



EE[{E - EY - ql] 
5|AtP tanh(/3£;/2) (2|AtP + - EE){E + E)[l - n{E) - n{E)] 



(51) 



(2|A, 



E{AE^ ~ ql) 
+ EE){E - E)[n{E) - n{E)] 



EE[{E + EY-ql] 



EE[{E - EY - ql] 



M{qo) = 0. 



(52) 
(53) 



It is straightforward to check that H{qo) — L{qo) ~ —1 and K{qo) ~ at large qo » |At| (here we dropped all real 
and imaginary contributions, suppressed by an inverse power of ql). In fact, these asymptotes are independent of 
the specific values of the temperature and the color superconducting gap, provided that the inequality T, \At\ ^ qo 
is satisfied. In this limit, therefore, the dispersion relations for both the magnetic and electric collective modes, see 
Eq. (E^) and (Ej), take the following form: 



for T, |At| < go < Z^- 



(54) 



where Up = g/i/\/27r plays the role of the plasma frequency. As we stated earlier, this result is in agreement with the 
assumptions used. 

For completeness of presentation, let us also write down the expressions for the component functions of the polar- 
ization tensor in limit of large go (and |g| =0). By making use of the definitions in Eq. (|2^), we derive 



Iliiqo) = U2{qo) u;] 
Il3{qo)^UM-0. 



(55) 
(56) 



The existence of the plasmon collective modes with large frequency is hardly a surprise in the dense quark matter. 
As in the case of ordinary metals, the plasmons are present in the normal and the color superconducting phases, and 
their properties are weakly affected by a small (compared to the chemical potential) temperature of the syst em an d 
by the presence of a gap. The "light" plasmons, that we reveal and study in this section too (see Subsec. VII E| ) , 
are much more interesting collective modes of the color- flavor locked phase of dense quark matter. In contrast to the 
ordinary plasmons, there is no trace of such "light" plasmon modes in the normal phase. Before studying them in 
detail, it makes sense, therefore, to review the properties of all collective modes in the normal phase first. 



A. Normal phase, T <^ ^ 



In the normal phase, the properties of the collective modes are well known. Here we just briefly review them. 

By definition, the value of the gap in the quasiparticle spectrum around the Fermi surface is zero in the normal 
phase. Thus, we substitute |Ay| = into our general expressions for the component functions in Eqs. (ES) - (pw. By 
further approximating the expressions in the limit of low temperatures (T <C /i) and small momenta (fgoli I'fl^ m)i 
we arrive at the following results: 

Hiq)^-l 



go go + |g| 
m 90 - k1 



(57) 
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K{q) 

m 



3 1 



90 , 90 

m 

90 



|91 



1 ^ ITTJ In 
2|91 90 



90 - \q\ 
90 + \q\ 



|91 
9o + I9I 



^ ^ i9i V m 90-I9I 

In the derivation, the ratio 90/I9I was kept arbitrary. Now, by making use of the definitions in Eq. (12 
well known hard-dense loop expression for the polarization tensor ||3q-|37|| : 



n,iq) = -Iltiq), 

n3(g) = n4(g) = 0, 



n;(9), 



given in the same representation as in Ref. 0|, with 



90 



90 + I9I 



'''^^^='"^1^^" 90-|9l 



(58) 
(59) 
(60) 

we derive the 

(61) 
(62) 
(63) 

(64) 
(65) 



In accordance with the approximation {T <^ /i), the finite temperature corrections were neglected here. Such correc- 
tions are suppressed by at least a factor of (T/^)'^. When the temperature is not smaU, one would need to correct 
the above result by also adding the contribution of the so-called hard thermal loops |3^ . 

It is straightforward to check that, in the limit \q\ <^ \qo\, this HDL polarization tensor indicates the existence of 



the plasmon collective modes with 



In the opposite limit, we reveal Debye screening (with mjj = S^^p) for the 



low frequency electric gluon modes, and no static screening for the magnetic gluon modes (there is, however, dynamic 
screening due to Landau damping). 



B. "Light" plasmon mode 



As we briefly mentioned earlier, the color flavor locked phase of dense quark matter reveals additional "light" 
plasmons that have no analogues in ordinary metals. In this subsection, we study these new modes at zero and finite 
temperature. We restrict ourselves to only the case of the long wave len gth limit. 

At r = 0, the expressions for the component functions in Eqs. ( ^l|) - (|53| ) further simplify: 



H{qo) = L{qo) 



- / de 

9 



9|Ao| 



K{qo) 
M{qo) 



E{4E^ - ql) 
, (2|Aop 



{2\Ao\^ - + EE){E + E) 
EE[{E + EY-ql] 



£2 - EE){E + E) 



E{AE^ 



2\ 
95) 



EE[{E + Ey 



9o^] 



= 0, 



(66) 

(67) 
(68) 



where |Ao| is the value of the color superconducting gap at zero temperature. 

The spectrum of massive (both electric and magnetic type) collective modes, then, is determined by the equation: 



ql + ujlHiqo) = ql+ ujlL{qo) = 0. 



(69) 



As is clear, there is at least one solution to this equation which corresponds to the ordinary plasmon mode, found 
earlier. We remind that such a solution appears at \qo\ ~ ^ |Ao|, where -ff(qo) = L{qo) ~ —1. Since the imaginary 
part of H(qo) = L^qo) is small at large |go|, this plasmon has a narrow width. 

It is more interesting to notice, however, that there is another mode with the energy less than a threshold of 



producing a pair of quasiparticles, |go| < 2|Ao|. Since |Ao| <^ Wp, the corresponding solution to Eq. (|69|) with a 
sufficient accuracy is given by a solution to the approximate equation, L(qo) — 0. Using Mathematica, we get the 
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solution, qo = tua ~ 1.362|Ao|. The approximate value of the mass of this mode could also be extracted from the 
derivative expansion of the effective action, obtained in Ref. [^. Using our approach, the corresponding derivative 
expansion appears after expanding the component function in powers of go: L{qo) « L{0) + (go/|Ao|)^ii(0) + . . .. It 
is straightforward to check that 



m 

Li(0) 



21-81og2 

54 ' 
21 + 16 log 2 

324 ' 



so that 



go = rriA 




21 + 16 log 2 



1.70|Ao| 



(70) 
(71) 

(72) 



We see that the derivative expansion somewhat overestimates the value of the mass. Still the result is below the 
threshold. Notice that value of the mass obtained in Ref. p3 was claimed to be go ~ 2.94|Ao| which is larger than 
our result by a factor of v^- We believe the authors of Ref. |39|] made a simple arithmetic mistake when substituting 
their expression (61) into Eq. (62). 

By calculating the component functions in Eqs. (^l|) - (|5^) numerically, we also obtained the temperature depen- 
dence of the light massive gluon mode. The result is presented in Fig. ^ 





FIG. 1. Left: the temperature dependence of the mass of the "light" plasmon measured in units of Tc . Right: the temperature 
dependence of the ratio mA/|AT|. 



VIII. NG BOSONS 



We recall that the original chiral symmetry is broken in the CFL phase of dense QCD. This means that an octet 
of pseudoscalar NG bosons should appear in the physical spectrum of our model. Following the approach used in 
Sec. we would like to study the axial current-current correlation function. To this end, we should introduce axial 
external sourses Af'^^ into the microscopic action. Then, after integrating out the quark degrees of freedom, we would 
get an effective action similar to that in Eq. (^3|), except that the quantum part of the gauge field would be absent. 
The subsequent integration over the NG boson field would lead us to a quadratic form similar to that in Eq. (|34|), 
determining the axial current-current correlation function: 



A -B 



u/q 



^5 



AB 



niOW(g) 



n2n3 + (n4 

Ha 



-0(l)(g) 



(73) 
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The location of the pole in this correlation function defines the dispersion relation of the NG bosons (which coinsides 



with the dispersion relation q^IV^^q'^ = 0, derived in Ref. |Q at zero temperature)]^ 

g^n^.q" = q^Hs = u^l [qlK(q) - \q\^ L{q) - 2goklM('?)] = 0. 



(74) 



Here we took into account the representation in Eq. ( |27q ). In the following two subsections, we analytically study 
this relation at zero temperature and in the nearcritical region. In the latter case, we also supplement our analysis 
by numerical calculations. 



A. Nambu-Goldstone modes, T = 

At first, let us consider the spectrum of the NG bosons at zero temperature. In this case, the component function 
M{q), see Eq. ( |3l] ) for the definition, is identically equal to zero. As is clear from Eq. ([74|), then, the inverse propagator 
of the NG bosons reads 

iV-^'^{q) = ^ [qlK{q) - |gf L(g)] . (75) 

In the far infrared region (Igojj I?! ^ I^oDj it is sufRcient to use the zero momentum limit of the functions K{q) and 
L{q). From Eqs. ( |66| ) and (|6^), we easily derive these quantities, 

ff(0) = m = ,76) 

= ,77) 

Finally, by making use of them, we obtain the well known expression for the decay constant as well as the low-energy 
dispersion relation of the NG bosons [^,^^ : 

^2 21 -8 In 2 , n 1,^2 

Png = -^ 18 ' '"''^ ^o-gkf- (78) 

The above dispersion relation remains almost unchanged even at small non-zero temperatures. In parti cular, in 
Appendix^, we show that the NG boson excitations acquire only an exponentially small width, see Eq. (EIS).^ With 
increasing the temperature, this width is expected to grow, reaching its maximum in the nearcritical region. 



B. NG bosons in the nearcritical region 



The existence of the small parameter \IS.t\/T in the nearcritical region is very helpful in applying analytical tools to 
studying the properties of collective modes. In the case of the NG bosons, whose spectrum is determined by Eq. (|74|), 
we need to calculate the component functions K{q), L{q) and M{q) as powers expa nsions in \/S.t \/T. To this end, it 
is convenient to utilize the representation of the component functions in Eqs. (50b) - Eqs. ( |50c|) . The approximate 
expressions for li integrals, with all the corrections up to second order in |At|/T, are given in Appendix ^ By 
making use of them, we derive the following relation: 



ng 

3 



-2 - 



57r|A2 
8T 



57r|A2 
ST 



1 - 



'ng 



2v, 



ng 



Vng . 1 

— - arcsm — 

2 t). 



ng 



49C(3)|A7 
47r2r2 



7C(3)|At| 
127r2T2 



(79) 



^We would like to thank V.A. Miransky for suggesting this alternative deivation of the NG boson dispersion relation. 

■^Notice that the dispersion relation of the NG bosons (and their width, in particular) might have corrections due to the 
interactions of the NG bosons with one another. Since their decay constant is of the same order as the chemical potential ^, 
such corrections are expected to be parametrically suppressed by an inverse power of ^. The detailed analysis of this point is 
beyond the scope of this paper. 



14 



where Vng = qo/\(l\ and C(3) « 1.202. The corrections to this equation are of third order in \ At\/T. Now, assuming 
that Vng vanishes as a power of |At| when T Tc, we arrive at a simple quadratic equation that determines the 
dispersion relation for the NG bosons, 



,2 7C(3)|AH^ _ 



.57r|A7 



(80) 



The solution to this dispersion equation reads 



qo = 



(81) 



where x'^g ~ 0.215 and j/*^ w 0.245. Notice that this solution is consistent with the assumption used in derivation of 

Eq. m. 



6-io' 





q=0.OOOG 
q=0 . 0010 
q=0 . 0014 



FIG. 2. Spectral density of the NG boson at T = 0.9997 (left) and at T = 0.9999 (right) as a function of energy. Everything 
is measured in units of Te- 



lia order to cross-check our analytical studies, we have also performed a numerical analysis of the NG boson spectral 
density. Such a density is defined by the imaginary part of the corresponding propagator. In the numerics, we have 
used the approximation of small energy and momenta, but we have not used any assumptions about the value of the 
superconducting gap. The results are plotted in Fig. || for two values of temperature, T — 0.9997rc (corresponding 
to |At| = 0.053Tc) and T = 0.9999Tc (corresponding to |At| = O.OSlTc). As is easy to check, the energy location of 
the maxima scale with momenta in accordance with a linear dispersion law that is very similar to that in Eq. (^l|). 



IX. CG MODES 



Now, let us consider a special type of collective modes, the so-called Carlson-Goldman gapless modes. Such modes 
were experimentally discovered by Carlson and Goldman about a three decades ago Q (see also, Ref. pl| , p0| -p5[ ) . 
One of the most intriguing interpretation connects such modes with a revival of the Nambu-Goldstone (NG) bosons 
in the superconducting phase where the Anderson-Higgs mechanism should commonly take place. It was argued 
that an interplay of two effects, screening and Landau damping, is crucial for the existence of the CG modes 
These modes can only appear in the vicinity of the critical temperature (in the broken phase) , where a large number 
of thermally excited quasiparticles leads to partial screening of the Coulomb interaction, and the Anderson-Higgs 
mechanism becomes inefficient. At the same time, the quasiparticles induce Landau damping which usually makes 
the CG modes overdamped in clean systems. To suppress such an effect and make the CG modes observable, one 
should consider dirty systems in which quasiparticle scatterings on impurities tend to reduce Landau damping. 

In the two fluid description, the CG modes are related to oscillations of the superfluid and the normal component 
in opposite directions . The local charge density remains zero in such oscillations, providing favorable conditions 
for gapless modes, in contrary to the widespread belief that the plasmons are the only collective modes in charged 
systems. 

The CG modes exist only in the presence of a large number of thermally excited quasiparticles (normal compo- 
nent). This means, therefore, that such modes could live only in a finite (possibly very small) vicinity of the critical 
temperature. 
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In the nearcritical region, from Eq. (^4[) we derive the following approximate equation for the spectrum of the 
collective modes of electric type (assuming that \/S.t\/T <C 1 and Vcg <C 1 where Vcg = Igol/M): 



2 457r|AT 



eg 



9C(3)|AtP 
87r2r2 



= 0. 



(82) 



The derivation of this relation is similar to the derivation of Eq. (|8C|) in the case of the NG bosons. It might look 
surprising, but Eq. (p^) is qualitatively the same as Eq. (p^). The solution to this dispersion equation reads 



90 



|At| 

T 



{±x* ~ iy* )\q\, 



(83) 



where x*^ « 0.193 and y*g sa 0.316. This solution corresponds to the gapless CG modes, and it closely resembles the 
dispersion relation of the NG bosons in Eq. (|8l|). The width of these CG modes is quite large, but this should not be 
surprising because we consider the clean limit of the dense quark matter. 

The numerical solution for Vcg as a function of temperature in Fig. |^. To plot the figure, we used the standard 
Bardeen-Cooper-Schrieffer (BCS) dependence of the value of the gap on temperature, obtained from the following 
implicit expression: 



In 



de[l- tanh ^ y?T]A7p 



(84) 



where 7 « 0.567 is the Euler constant. As was shown in Ref. |Tj], such a dependence remains adequate in the case of 
a color superconductor. 




0.988 



0.992 



0.996 



T/Tn 




FIG. 3. Left: The real (solid line) and imaginary (dash line) parts of the numerical solution for v = qo/\Q\- Right: Temper- 
ature dependence of the color superconducting gap. 



It has to be emphasized that the result in Eq. ( |83| ) is obtained for a clean system where Landau damping effects 
have their full strength. Therefore, it should have been expected that the CG modes are overdamped. However, our 
analysis shows that the presence of two different types of quarks with nonequal gaps in the CFL phase plays the key 
role in a partial suppression of Landau damping. Indeed, as is straightforward to check, if the values of the gaps in 
both the octet and singlet channels were equal, the dispersion relation of the CG modes would be given by 



. 14C(3)|At| , 
37r3T ^ 



(85) 



where the higher order terms in powers of \/^t\IT are denoted by the ellipsis. We see from this result that the ratio of 
the imaginary and real parts is much larger than 1 (in fact, it is infinite in this leading order approximation), meaning 
that the CG modes are unobservable in the case of equal gaps. In contrast to this, the CG modes with the dispersion 
relation in Eq. (p^) could be observable (see our results for the electric gluon spectral density below). 

From the numerical result, we clearly see that the ratio of the width, F = —2Im{vcg)\q\, to the energy of the CG 
mode, Eq = Re{vcg)\q\^ increases when the temperature of the system goes further away from the critical point. At 
temperature T* « 0.986Tc, this ratio formally goes to infinity. This value of T* gives an estimate of temperature 
where the CG mode disappears. Of course, such an estimate is not very reliable since our approximations should 
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break before the temperature T* is reached. In order to get a better estimate of temperature where the CG mode 
disappears, a more careful numerical calculation is required. 




As in the case of the NG bosons, we have also performed a numerical analysis of the electric gluon spectral density 
in the region of small momenta where the trace of the CG modes is expected to appear. In our numerical calculation, 
we have used the approximation of small energy and momenta. The results are plotted in Fig. ^ for two values of 
temperature, T = 0.9997rc (corresponding to |At| = 0.053Tc) and T = 0.9999Tc (corresponding to |At| = O.OSlTc). 
Although the maxima in the electron gluon spectral densities are not so well pronounced as in the case of the spectral 
densities of the NG bosons (see Fig. ^ in the preceding section) , they are detectable when temperatures are sufficiently 
close to Tc- One could also check that the energy location of the maxima in Fig. |4| scale with momenta in accordance 
with a linear dispersion law similar to that in Eq. (p3|). 



X. MEISSNER EFFECT 

Now, let us discuss the Meissner effect at finite temperature. This can be done by examining the response of the 
quark system to an external static magnetic field. In order to derive the expression for the magnetic ^i{q) component 
of the polarization tensor in the limit jgol ^ I*?! — * 0, let us start with the calculation of the iJ(go — 0, function. 



diil - n / de 



2\An 



-1 

.2 I 9'?^ 



2|A. 



7e^ 



79_£_ 
' 4 



2|A. 



3l£ 



El) 



{u:l + El){u:l + EX) 



{.^l + El){i^l + El) 



(86) 



where E\ = (e± q^/2)^ + I^Tp and E\ = (e± q^/2)^ + 4|AtP. The proper calculation of this expression requires 
performing the Matsubara summation before integrating over e. This order could be interchanged if the normal phase 
contribution (with Ay = 0) is subtracted from the above expression, i.e.. 



T 



Hm = H,{\q\) + [H{\q\) ~ H,m] = ffo(ki) + E / ^^(1 - ^' 



de 



2|A. 



7e^ + 7^ 



{u^l + Ei){u^l+Ei) 



2IA7 



4 



2IA. 



4 



{ojl + El){u:l + El) 



{ul + El){u:l + E^, 



- 9 



4 



{u:l + {e-q^l2Y){u:l + {e + q^/2Y) 



(87) 



It is easy to check that performing summation over Matsubara frequencies (and taking the limit |q| 0) we come 
back to Eq. (^6|) (where we put first go = first and then |(?| = 0}. Notice that the last term in square brackets in 
Eq. ( p7| ) reproduces the first term on the right hand side in Eq. (^6|). 

The nice feature about Eq. ( |8^ ) is that it allows to interchange the order of summation and integration. This is 
because the integral and the sum converge rapidly. The initial expression for if(qo, |g|) did not have this property, and 
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the order of operations was fixed (first summation, then integration). By evaluating the integral over e in Eq. (|8 
we obtain 



ttTIAtI 
12~~ 



5/2 



9|At|4 + 2(72^2(2^2 + 5|Aj,|2 + ^2^2/3) 



+ |AT|2 + g2^2/4) 

- 3|AtP + 9^^2/3 2(c.2 + 3|AtP + q^e/^) 



V^?. + |At| 



v/^^ + 4|At|2 



From the structure of this expression, it is clear that |(z| /Tc is an important dimensionless parameter at all temper- 
atures. Indeed, for T <C Tc, it is the gap |Ao| which of order Tc that sets the scale in the H{\q]) function. Near Tc 
(i.e., \Tc — T| <C Tc), on the other hand, the gap is small but ujn^o — 2miT ~ Tc gives the same characteristic scale. 

Let us begin by examining the first limiting case, \q\ <C Tc- In this limit, we keep only the leading order term on 
the right hand side of Eq. (pq), and drop all terms suppressed by a power of Thus, we arrive at 



H{\q\) 



ttTIAtI' 



-5/2 



-3|At|V2 



(C^2 + |At|2)3/2 3|A^|4 I ^^2+|A^|2 ^ 



(89) 



As a benchmark test, let us calculate this expression at T = 0. The sum over frequencies in this case can be replaced 
by an integral, by setting 2T:TSn = duj. Thus, we arrive at the result 



HM - 0) 



lAi 
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duj 



-5/2 



LO 



-3|At|V2 



3|AtP 



dx 



(a;2, 
5/2 



(2.2 ^ 1)3/2 



|At|2)3/2 



3|At|4 \^ ^c^2 + |Aj,|2 ^:j2q 

3/2 a;2_^3 \ i 21-81og2 



+ 1 Va;2 + 4 



54 



(90) 



which coincides with the known expression at zero temperature [£2|,^,^6|3J] . 

Near the critical temperature Tc, we expand Eq. (89) in powers of |At| to obtain 



27rT|AT 



2 °o 



E „,3 



ji=0 



7C(3)|At| 
127r2T2 



(91) 



The nonzero value of i?(|(7| — > 0) function in the static limit indicates that a constant magnetic field is expelled from 
the bulk of a superconductor. So, the conventional Meissner effect takes places at all temperatures T < Tc in the 
color superconducting phase. 

For completeness of our presentation, let us also determine the penetration depth of the magnetic field. To this 
end, we also need to calculate the behavior of the i?(|(?|) function for large values of Thus, let us consider the limit 
\q\ ^ Tc. In this case, the main contribution to the integral over ^ in Eq. ( p8| ) comes from the region ^ < Tc/\q\ ^ 1. 
In view of this observation, it is justified to neglect —^2 term in the overall factor (1 — ^2) Then, by substituting 
[(fl^ = a; and expanding the limits of integration to infinity, we obtain the following approximate result: 



7r2T|AT|' 

3|g1 



E 



5/2 



1 



|AtP 
3 



|At|2v/c.2+4|a^|2 



or, written in slightly different form, 



HM) 



-1 



4|A2 



■ tanh ■ 



;2 -K4|ArP + \/w2 + |Ay 



^E 



F) 



(92) 



2T 



9|A2 



^„ V^2 + |At|2v/c.2+4|A^|: 



Vc.2+4|Ar|2 + ^C.2-t-|AT|2 



(93) 
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Here we explicitly performed the summation in first term and slightly rearranged the second term. In calculation 
leading to Eq. (|9^), we used the following table integrals: 



dx 



— oo — c 

At zero temperature the expression (|9^) reduces to 



dxx^ 



2/3x2 + 1 V2\/^TT' 



(94) 



H{\q\) 



TrjArl 
6|91 

ttIAt 



T 



dx{2 + x^ - \Jx^ + 4Va;^ + 1) 







Va;2 + 4\/a;2 + 1 



Near Tc, we expand the expression in Eq. (|93| ) in powers of |At|, and obtain 



(95) 



(96) 



Now, we consider the calculation of the magnetic penetration length. By making use of the standard definition (see, 
for example, Ref. |p6[), we obtain the following estimate: 



^=1 

TT 



dq 



g2 - ujlHiq) TT 



l/?o 



dq 



g2 + m\f 



dq 



l/Co 



q^ + M^\AT\/q 



(97) 



where we split the range of integration into two qualitatively different regions with momenta smaller and larger than 
the coherence length ~ l/^c, respectively. The parameters mj^j and M'^ are defined as follows: 



m 



M 



^pT^' for T«r„ 



7C(3)|ATr ,, ,2 f„„ T T 

Mo2 = f + 2^(V3/2) - if (V3/2)) ^2, for T « T„ 

for T ~ Tp. 



^|At| , ,2 
8T "^P 



(98) 
(99) 



After performing the integration in Eq. (97), we derive 
6 = 



"^M + 575(^^'|AT|)-l/^ for (M2|AT|)i/3»l/eo, 
mil, for (M2|Aj,|)i/3«i/^„. 



(100) 



At low temperatures, T <^ Tc, we have 1/Co ~ |At| and ~ ^^o^ ^ '^i where (5l is the London penetration 
depth. By taking into account that Mq LUp ^ |At|, we arrive at the following low temperature expression for the 
penetration depth of a magnetic field: 



S = Sp = 



3V3 



(A^2|AT|)-l/^ r«r. 



(101) 



This is the so-called Pippard e xpre ssion for the penetration depth (note that i5p ^ ^o)- It is easy to check from 
the general expression in Eq. (IOC) that S = 5p almost in the whole region of temperatures T < Tc- It is also 
straightforward to check that the London limit (by definition, 5 ^ ^o) is realized only in a very close vicinity of 

the critical line where \Tc-T\ < Tc. 
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XI. CONCLUSIONS 



In this paper, we studied collective modes coupled to (vector and axial) color current in the color-flavor locked 
phase of cold dense QCD at zero and finite temperatures. In this class of collective modes, we revealed and studied 
the ordinary plasmons, the new type of "light" plasmons, the NG bosons and, finally, the gapless CG modes. 

As we argued in the main text, the properties of the ordinary (high frequency) plasmon collective excitations are 
similar to those in ordinary metals. The plasma frequency LUp is proportional to the value of the chemical potential, 
and it is essentially independent of the values of the temperature and/or the superconducting gap in the dense quark 
matter. The "light" plasmon, on the other hand, has very different properties. It is also a massive excitation (that is 
why we also call it a plasmon), but the value of its mass ttia is of the order of the superconducting gap (more precisely, 
1.362 1 At I < ruA < 2|At|) which is quite small compared to ujp. The appearance of such a mode seems to be directly 
related to the existence of two different gaps in the quasiparticle spectra (in the octet and singlet channels). 

In the long wave length limit \q\ — > 0, the "light" plasmons are stable with respect to the decays into the quark 
type quasiparticles. It is only because of other possible decay channels (for example, involving NG bosons), that 
their width might be nonzero. Now, while we have not studied the detailed dispersion relation of the "light" plasmon 
modes in the short wave length limit \q\ ^ IAt'I, it is expected that their energy is a monotonically increasing function 
of momentum \q\. If this is true, there should exist a critical value of the wave length of order 1/|At| at which the 
energy of the "light" plasmon becomes equal to the threshold of the quasiparticle pair production. The excitations 
with the wave lengths shorter than the critical value (and the energy larger than 2|At|) could easily decay into 
quasiparticles and, therefore, should have a relatively large width. It means that these new type plasmon modes have 
two characteristic plasma frequencies, ttia and 2|At|. The stable (narrow width) modes live only in the following 
window of the energies, ttia < qo < 2|At|. 

The properties of the NG bosons in the CFL phase at zero temperature are well studied in the literature. Here we 
generalized such studies to the case of finite temperatures. In particular, we showed that the NG boson properties 
at small temperatures remain almost the same as at zero temperature, having an exponentially small width. With 
increasing the temperature, the width also increases. We presented analytical study of the dispersion relation of NG 
bosons in the nearcritical region of temperatures. Our result confirms the general expectation of the slowdown of the 
NG bosons in this region, where their maximum velocity is proportional to the small value of the color superconducting 
gap. 

By making use of an explicitly gauge covariant approach, we also studied the properties of the gapless CG modes 
in the CFL phase of cold dense quark matter in the nearcritical region (just below Tc) where a considerable density 
of thermally excited quasiparticles is present. It is important to mention that the presence of the CG modes coexists 
with the usual Meissner effect, i.e., an external static magnetic field is expelled from the bulk of a superconductor. The 
existence of the Meissner effect is a clear signature that the system remains in the symmetry broken (superconducting) 
phase. 

In the case of the CFL phase, as we showed, the CG modes appear even in the clean limit. Despite the sizable 
width, their traces can be observed in the spectral density of the electric gluons. Taking the effect of impurities into 
account should, in general, make the CG modes more pronounced [Q. In realistic systems such as compact stars, 
natural impurities of different nature could further improve the quality of the gapless CG modes. 

The existence of a gapless scalar CG modes (in addition to the pseudoscalar NG bosons) is a very important property 
of the color superconducting phase. They may affect thermodynamical as well as transport properties of the system 
in the nearcritical region. In its turn, this might have a profound effect on the evolution of forming compact stars. 
The CG modes might also have a consequence on a possible existence of the hypothetical quark-hadron continuity, 
suggested in Ref. |47|. Indeed, one should notice that, in the hadron phase, there does not seem to exist any low 
energy excitations with the quantum numbers matching those of the CG modes. 

In the future, it would be interesting to generalize our analysis to the so-called S2C phase of dense QCD. In 
absence of true NG bosons in the S2C phase, the gapless CG modes might play a more important role. Our general 
observations suggest that Landau damping should have stronger influence in the case of two flavors. At the same 
time, the presence of massless quarks could lead to widening the range of temperatures where the CG modes exist. 
To make a more specific prediction, one should study the problem in detail. 
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APPENDIX A: IDENTITIES, FORMULAS, ETC. 

In the main text we use the following identities for the structure constants of SU{3): 

^ACD^BCD ^ ^ ,AB^ (^1) 
o 

fACDfBCD ^^^AB (^2) 

For completeness, we also present the following more complicated identities: 

(A3) 
(A4) 
(A5) 
(A6) 

In this Appendix, we also derive the three types of Matsubara summation formulas that appear in the main text. 
By convention, we use the notation a;„ = T7r(2n+ 1), = 2T-Km and /3 = 1/T. The mentioned three types of sums 
are 

"'■''""^''■'''°^? K + »')[K + o,„)^ + i,T 

a, J-) = + (A9) 

Let us start from the first one. The sum is performed by going to the complex plane and introducing the contour 
integral. The result is 

_ 6(0^ - + &') tanh(a/?/2) + a(0^ + ~ b^) tanh(6/3/2) 
a, b) - [{a + bf + [{a bf + 

_ (a + b) [1 - 11(a) - ii{b)] (a - b) [n{a) - n{b)] 
2ab [{a + 6)2 + ^1] + 2ab [{a - bf + ' 

where 

n{a) = — - (All) 

^ exp(a/3) + 1 ^ ' 

is the Fermi distribution function. 

To perform the analytical continuation to the real frequencies, we substitute iClm = O + iry where is a vanishingly 
small positive constant. In this way, we arrive at 

ri[ii,a,0)- y.V. 2a6[r!2_(a + 6)2] ' ' 2ab[Q? - [a - bf] 
+z^^^|-t^sgn(n)(5 (1^2 _ (a + 6)^) [1 - n(a) - n{b)] 

+i '^^°' \^\ gn{Q.)5 (02 _ (a _ hf) [n{a) - n{b)] . (A12) 



(AlO) 



21 



Similarly, 



iQm [2a6tanli(a/?/2) - (a^ + b^ + O^) tanh(6/3/2)] 

n 

F3{i^m;a,b) 



2b [{a + by + f]2j 2b [{a - by + ftl] ' 
a(a2 - b- + Vl,){-Mi\i{(i.l/2) + h{lr - ir + 0;J t;uili(/AJ/2) 

2[{a + by + nu [{a-hy + nu 

(g + b) [1 - n{a) - n{b)] (a - b) \n{a) - n(b)] 
2[(a + 6)2 + n^] 2 [(a -6)2 + 0^] ' 



and 



.tt\^\ 



2b 

F,{Q, a, b) = -V.V. 2[J72-(a + 6)2] + 2 [f^^ _ (« _ 6)^] 

.Tr{a-b) ,r,-,r/r,2 



2 

The typical expressions for a and b would be given by 



where e = |p| — /U, ^ = cos 6 and 9 is the angle between the spatial vectors p and q. 

The (5-functions that appear in the imaginary terms of all three sums have the singularities at 



,± 





|g]^(Ai - Al)±n^[n^ - - (Ai - A^yW - we - (Ai + A^FT 

2(02 _ 1^-12^2) 



(A13) 



(A14) 



26 [Q2 _ (a + 6)2] ■ '26 [02 -(a -6)2] 
i^S _ (a + 6)2) [1 _ n{a) - n{b)] 

S (0^ - (a - bf) [n{a) - n{b)] (A15) 



sgn(0)5 (0^ - (a - 6)^) [n(a) - n(6)] . (A16) 



a = ^A2 + (e-e|g1/2)2, (A17) 



6= jA2 + (e + ^1^1/2)2, (A18) 



(A19) 



which is real only when 02 < (Ai — A2)^ + or when 02 > (Ai + A2)^ + I^I^C^- The first region defines the 
window of energies in which the Landau damping operates, while the other region corresponds to energies for which 
the on-shell pair production of quasiparticles is allowed. 

APPENDIX B: OPERATORS 0<^> 

In the main text of the paper, we heavily use the following set of four tensors: 

Ol^,^ (q) = 9tJ.iy - Uf,u^ + ^j^, (Bl) 

0'-^M=u,u.-'^-'^, (B2) 

O^^Hl) = '-f^, (B3) 

0^,tHa)=02.^| + |.Mt^. (B4) 
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The first three of them are the same projectors of the magnetic, electric and unphysical (longitudinal in a 3+1 
dimensional sense) modes of gluons which were used in Ref . ||7| . In addition, here we also introduced the intervening 
operator O^t^(g), mixing the electric and unphysical modes. In the above definitions of 0^^l{q), we use 



u^ = (1,0,0,0), 

By making use of the explicit form of the operators O^'' , it is straightforward to derive the following relations: 



^^--^o(2.H,) + |o(l)(g) 



|g1 q' 

as well as the following multiplication rules: 



o^'Jiq)~o^^Xq) 



qo\q] 



ql 



(B5) 
(B6) 



(B7) 
(B8) 



(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
(B17) 
(B18) 

In these last expressions, proper contractions of the Lorentz indices are tacitly assumed. 

The structure of the polarization tensor as well as the structure of the gluon propagator in the main text are given 
in terms of the O*-*-* operators. It is ver y u sefu l, therefore, to know how to invert the corresponding tensors. By 
making use of the definitions in Eqs. (Bl) - (B4), we check that the following statement is true. If a general matrix 
(in Lorentz indices) allows the following decomposition in terms of the operators O'^^^ 
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(g)OW 
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2)((7)-0(3) 
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(q) 
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(q) 




{q)0^^\q). 




0(4) 


(q)0(3) 


(q) 


= 0(2 


{q)0^^\q). 





A = aO(i) (g) + 60(2) (g) + cO^^) + dO^^) (g) 

its inverse matrix is given by 

A-i = ioW(g) 



bc + cF 



0^'\q) 



bc + (P 



bc + d^ 



O^'Hq). 



(B19) 



(B20) 



APPENDIX C: CALCULATION OF U^^ 

In the Appendix, we calculate the expression of II^j^ in different limits. The complete expression for the polarization 
tensor contains a sum of the left- and right-handed contributions [only one of them is given in Eq. (|l|) in the mam 
text] : 

K^{q) = ^tUi) + ^iUi) ^ S^^'^.Aq)- (CI) 

The complete one-loop result (in Minkowski space) is given by 
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|ee'ep_q/2<^p+9/2'^M^'' 7Fi(ri; £'p_,/2, -E'p+g/2) + ^li^'^ ^ 
-T^l"" 7F3{Q; -Bp_g/2' + ^si^: ^ 



'p-g/2' ^p+9/2) + E^i^; £^p-g/2i ^p+9/2) + 



+2|At|V;/ 



_,/2) +^2(^^;£^p_q/2,£^p+q/2 
7Fl(il; £'p_q/2, ■E'p+g/2) + ■E'p-9/2' ■E^p+g/2) + ■F'llf^; E'p~q/2^ E'p+q/2) 

Fli^i Ep_q/2^ Ep+q/2) + F^i^i Ep-q/2^ ^p+q/2) + Fli^i ^p-q/2^ Ep+q/2) ^ (C2) 



where the tensor structures are 



V-'' = tr 



(7M7oAp^g/27^7oAp_^|j/2) - (7m7oA;_,/27<^7oA;'+,/2) , 
' {^t^K-q/2^'^K+q/2) + tr (7/.A;_^/27.A^/^/2) ' 



and 



(e|)2 + |AT| 



(C3) 
(C4) 
(C5) 

(C6) 
(C7) 



are the excitation energies in the oct et and the singlet channels, respectively. 

The general expression in Eq. (C2) can be approximated by dropping all the subleading terms. For large chemical 
potential, we derive the following approximate expression: 

Uplift, q) = {Uf,u„ - g^,^) 



Fi{n;E^,E+)+Fi{n;E^,E+) + Fi{n;E^,E+) 



2 "f"" 

7Fi {Q; E^,E+)+Fi (^i; E^,E+) + Fi (^i; E^,E+) 
^e^M^ + u^q.) [7F,in; E_,E+) + F,{n; E_,E+) + F,{n; E_,E+) 
u.u^l + e) eg,. ^^OW(g) 
7F3{Q;E^,E+)+F3in;E^,E+)+F3in;E^,E+) 



2eC 



k1 



7F2(f7; E-,E+) + F2in; E- , E+) + F2(f7; E- , E^ 



(C8) 



APPENDIX D: KEY INTEGRALS IN THE NEARCRITICAL REGION 



In this Appendix, we consider the nearcritical asymptotics of all the integrals entering the definition of the compo- 
nent functions H{q), K{q), L{q) and M{q), see Eqs ( |50a[ ) - (50d). It is sufficient for our purposes to consider only 
the Hmit of small momenta, IgoU?] ^ I^tI, keeping the ratio v — \qo\/\q\ arbitrary. By expanding the results in 
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powers of the ratio | Ayj/T and keeping the leading order terms (up to the second power in | A^l/T), we 
following expressions: 
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where K{z) and E{z) are the complete elhptic integrals of the first and second kind, C(z) is the Riemann zeta function, 
and F(a, 6; c; z) is the hypergeometric function. Notice that 
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APPENDIX E: THE RESULTS FOR h THROUGH /g AROUND T = 



In this Appendix, we consider the asymptotics of the integrals (see Appendix |d| for their definitions) at small 
temperatures. The approximate expressions for the first five integrals are straightforward to derive, provided T <C 
|At| « |Ao|. The resuhs read 
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The calculation of the other four integrals is more complicated. For example, we obtain the following representation 
for the /8(f) function: 
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Only the imaginary part of Is,{v) can be calculated exactly 
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where Ei{~z) is the integral exponential function. Assuming that 

1 



and 



T T VVT^ 

we derive the following asymptotic behavior of the imaginary part: 

Im/8(u) ~ e - 
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Similarly, we calculate the asymptotics for the functions Im/g(w) and Im/7(w): 



\mlQ{v) 
Im/7(ti) ~ — 2l7rti 



'pL 

V 



7-^(l-w2)(l + v/l-i'2) 



(E8) 
(E9) 

(ElO) 
(Ell) 



Neglecting for a moment exponentially small temperature corrections to the real parts, we can write down the 
dispersion relation for the NG bosons. 
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This gives the following solution: 
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Note that the dispersion relation of the NG bosons has an exponentially small imaginary part near T = (compare 
with the corresponding dispersion relation for the Anderson-Bogolyubov mode in ordinary superconductors |}42| ). 
Clearly, this qualitative feature of the result would not change even after (exponentially) small real contributions are 
added to h{v), l7{v) and /8(f)- 
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